High energy asymptotics and trace formulae for the 
perturbed harmonic oscillator 

■ Alexander Pushnitski* and Ian Sorrell^ 
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^ Abstract 

A one-dimensional quantum harmonic oscillator perturbed by a smooth compactly 
supported potential is considered. For the corresponding eigenvalues A n , a complete 
asymptotic expansion for large n is obtained, and the coefficients of this expansion are 
Oh ■ expressed in terms of the heat invariants. A sequence of trace formulas is obtained, 

expressing regularised sums of integer powers of eigenvalues A n in terms of the heat 
_C invariants. 
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1 Introduction and main results 

1. Local heat invariants. In order to state our main results, we need to recall the notion 
of local heat invariants. Let v G C°°(R) be any real valued function such that v and all 
^ | derivatives of v are uniformly bounded on R. For the self-adjoint operator h = — + v in 
L 2 (R, dx), consider its heat kernel e~ th (x,y), t > 0, x, y G R, i.e. the integral kernel of the 
■ operator e~ th . For any x E R, one has the asymptotic expansion 

-t— > 

Si » 
e- th { x ,x)^-—^a j [v{x)i t^+0, (1.1) 

X 

where aj[v (x)] are polynomials in v and derivatives of v, known as the local heat invariants 



of h. Explicit formula for aj[v (x)] is available: 



a A< x )\ = E m^^wkzT^ + <yn - y? k ) u • (i-2) 
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Formula (|1.2j) was derived in [7] on the basis of the results of Ell ; see also references in 
to earlier works on this subject. From (J1.2)) or otherwise, one obtains 



a [v(x)) 

a 3 [v(x)] 
a 4 [v(x)] 



1, ai[v(x)) 



v(x), a 2 [v(x)] = -v 2 (x) 
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2. Perturbed harmonic oscillator. Consider the self-adjoint operators 



o — 



dx 2 



+ x 2 and H 



dx 2 



+ x 2 + q(x) in L 2 (R, dx), where g G C£°(R). 



These operators can be defined as the closures of the symmetric operators, defined on 
Co°(R) by the same differential expressions. Denote by A° = 2n — 1, n — 1,2,... the 
eigenvalues of H and by Ai < A 2 < • • ■ the eigenvalues of H. The aims of this paper are 
(i) to describe the asymptotic expansion of A n as n —>■ oo, including explicit formulas for 
the coefficients of this expansion in terms of the local heat invariants; (ii) to derive a series 
of identities (trace formulas) which relate regularized sums of the type Y^m=i ^n, J2^=i 
etc. to some explicit integrals involving heat invariants. 

Our results are modelled on the Gel'fand-Levitan-Dikh trace formulae for the Sturm- 
Liouville operator (see 001313 or jlj) and in part motivated by the recent advances in 
calculation of the heat invariants [HI E3 • 

First, as a preliminary result, we establish the asymptotic expansion 



where cij are the local heat invariants. In formula ()1.3j) (as elsewhere in this paper) q G 
Cq°(R) and thus all the integrals in the r.h.s. converge. On the formal level, (jl.Hjl follows 
by subtracting (jl.lj) with v (x) = x 2 from (jl.lj) with v (x) = x 2 + q(x) and integrating over 
x. A rigourous justification of this formal procedure is not difficult and will be given in 
Section El 

3. High energy asymptotics. Suppose that q is given. Due to the explicit formula 
(jl.2j) . we can regard the integrals appearing in the r.h.s. of as known quantities. 

Below we describe the asymptotics of eigenvalues A n in terms of these integrals. Here is 
our main result: 

Theorem 1.1. (i) One has the asymptotic expansion 



Tr(e 



-tH 



1 °° f 

) =Vf J / (aj[x 2 + q(x)} - aj[x 2 ])dx, t +0, (1.3) 

V 4vrt Jr 




n 



oo, 



(1.4) 



with some coefficients Cj G R. 



2 



(ii) The coefficients Cj in (jl.4j) can be calculated in the following way. Consider the 
formal asymptotic expansion 

oo , 

^~^ + E7T%T' n ^°°' (L5) 



J- 



~1 i^n) 



with the coefficients 

bj = (V^r(| - j))- 1 [ ( aj [x 2 + q{x)\ - a^x^dx. (1.6) 
Jr 

Then inverting the formal asymptotic series (ll.5j) gives (|1.4)1 . 

Remark 1.2. 1. Theorem 11.11 gives an algorithm of computing the 'unknown' coefficients 
Cj in the expansion (jl.4j) in terms of the 'known' integrals (jl.6|) . The algorithm is given in 
the form of inverting an asymptotic series, which is a well defined algebraic procedure. 

In order to compute a coefficient Cj, one needs to know finitely many coefficients bj. For 
example, 

ci = c 2 = 0, c 3 = -b 2 , c 4 = --ft?, c 5 = -6 3 , c 6 = -2bib 2 , c 7 = -b\ - 6 4 . 

2 o 

2. The fact that only half-integer (and not whole integer) negative powers of A n are 
present in the r.h.s. of (jl.5j) is equivalent to a series of identities for the coefficients Cj. For 
example, the first three identities of this type are 

c 2 = 0, c\ + 2c 4 = 0, c 6 + c\ + 2cic 3 = 0. 



3. From Theorem II. lf ii) we obtain, in particular, 



1 r , , , 1 y . . 2 , 1 /" 2 

ci = — / q(x)dx, C2 = 0, C3 = — / aa; H / g 

Jr ^ Jr 27T y 



3/ J doc . 



1 1 Z" 00 1 

c 4 = --Cj, c 5 = — / (g 3 (x) + 3g 2 (x)x 2 + 3g(x)x 4 + -{q\x)f + 2q(x))dx. 

2 167T J.oo 2 

4. Trace formulas. As a by-product of our construction, we also obtain trace formulas 
for the eigenvalues \ n and A°. This result is a direct analogue of the trace formulas for the 
Sturm-Liouville problem due to 121 El IS] and our proof follows the reasoning of j2j • Let 
us introduce the Zeta functions 

00 00 

Z(s) = Y,K s , Z (s) = J2(^r s , Res>l. (1.7) 

n=l n=l 

If X n < for some n, then A~ s should be understood as |A n |~ s e~""\ If A„ = for some n, 
then the corresponding term in the sum Y2^Li ^n. s i s omitted. 

Due to the explicit formula \° n = 2n — 1, we have Z (s) = (1 — 2~ s )((s), where £(s) is the 
Riemann Zeta function. By the properties of (, we conclude that Z (s) has a meromorphic 
continuation into the whole complex plane with the only pole at s — 1, and this pole is 
simple. The real zeros of Zq are at s = —2n, n = 0, 1, 2, ... . 
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Theorem 1.3. The function Z(s) admits meromorphic continuation into the whole complex 

1 _3 _5 

2' 2' 2> 



plane. Its poles are simple and located at s = 1 and at s = — \, — \, — %, ■■■■ We have the 



identities: 

Z(-k) = Zo(-Jfe), keN. (1.8) 

As in |2 E] , formula (jl.8|) can be combined with the asymptotic expansion (jl.4j) to obtain 
regularised trace identities as follows. Exponentiating the asymptotics (jl.4|) . we obtain for 
any Res > 1: 

CO 

3=0 

where dj(s) are explicit polynomials in s and Cj. For example, 

d (s) = 1, di(s) = d 2 (s) = 0, d 3 (s) = —sci, d 4 (s) = -sc 2 , d 5 (s) = —sc 3 , 

d 6 (s) = -sc 4 + g ( g + ^ c\, d 7 (s) = -sc 5 + s(s + \)c x c 2 . 
Using this notation, we have for any fceN: 

oo 2fc+2 2fc+2 

n=0 j=0 j=0 

Now both sides of (jl.lOj) can be meromorphically continued into the half-plane Res > 
—k — ~. By Theorem II .3| the l.h.s. of (jl.lOj) is analytic at s = —k. By (jl.9j) . the same 
applies to the first term in the r.h.s. of (jl.lOj) . Thus, the second term in the r.h.s. of (jl.lOj) 
is also analytic at s = —k. As Z (s) has a pole at s = 1 (and no other poles), it follows 
that d 2 k+2(—k) = 0. Thus, we obtain 

oo 2k+l 2k+l 

Z(~k) = ~ E d i(- k )( X n) k ~ m } + E dA-k)Z {-k + (j/2)). 

n=0 j=0 j=0 

Combined with (jl.8j) . this yields a series of formulas 

oo 2fc+l 2k+l 

E< A « " E d i(- k K X n) k ~ U/2) } + E dA-^Zai-k + (j/2)) =0, keK (1.11) 

n=0 j=0 i=l 

In particular, for k — 1, 2, 3 we obtain (taking into account that Zq(0) = 0) 

oo 

E( A « - X °n ~ ^0=) + c i^o(-|) = 0; (1.12) 



n=0 



2r 

E( A ' - ( A °) 2 - 2c i V 7 ^ - -m) + 2 Cl ^ (-|) + 2c 3 Z (i) = 0; (1.13) 

n=0 V *n 

oo 

E K " ( A n) 3 " 3C l( A n) 3/2 " ~ 3(c 4 + C?) - 3c 5 (A°)-^) 

n=0 

+3 Cl Z (-|) + 3c 3 Z (-i)+3c 5 Z (i) = 0. 
Formulas (jl.l2jl . (j!.13j) (in a slightly different form) were obtained earlier in j^j. 
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2 Proof of Theorem ll.l( ii) 

The proof of part (i) of Tfieorem ll.ll is fairly standard and is based on the asymptotic theory 
of solutions to ODEs and on various explicit formulas for parabolic cylinder functions (which 
give the solutions to the ODE corresponding to q — 0). We give this proof in Sections 5-6. 
The proof of part (ii) of Theorem 11.11 is the core of our construction and is presented in 
this section. The proof is based on the following 

Lemma 2.1. Let A° = 2n — 1, n eN, and let Ai < A2 < • ■ ■ be a sequence of real numbers 
such that X n = X° n + 0(1) as n — ► 00. Suppose that an asymptotic expansion 

00 00 

3=1 3=1 

holds true, where ^ cti < 012 < ■ • • are some non-integer exponents and {pj} C M., 
{qj} C M. Then one has the asymptotic expansion 

E e " A " ~ Yt + E l r(1 - ^ + 2 logt £ ^T^rk^ + E ^ ( 2 - 2 ) 

n=l j j=l ^ k=l 

as t — > +0, wif/i some coefficients {r^} C R. 

Proof of Theorem 1.1(H). Given Lemma T2. II and part (i) of Theorem ll.il the proof of The- 
orem ll.lf ii) is immediate. Indeed, inverting the asymptotic expansion (jl.4j) yields the 
expansion of the form 

00 00 
\° n ~ A n + 2J 3 + X^ftjA" 3 ', n — > 00 

3=1 3=1 

with some real coefficients {bj}, {bj}. Now using Lemma 12.11 and the explicit formula 
Yl™=i e~ tx " = (2sinht) _1 , we obtain the asymptotic expansion 

n=l Vt i=l j=l U J ' fc=l 

with some real coefficients {jk}- Comparing this to (jl.3|) . we see that all coefficients bi 
vanish and the coefficients bj are related to the heat invariants by formulas (jl.6|) . This 
completes the proof of Theorem ll.lf ii). ■ 

In the rest of this section, we prove Lemma f2. II Broadly speaking, this Lemma can be 
regarded as a discrete analogue of the following version of Watson's Lemma: 
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Lemma 2.2. Let ip : R — > R be a locally bounded measurable function, such that ?/>(A) = 
for all A near — oo. Suppose that ip has the following asymptotic expansion 

where {aj} Cl\N, {f3j} C N, {pj} C R, {qj} C R are finite sets and M > max({a i } U 
{/3j}), M G (0, oo) \ N. Then the following asymptotic formula for the Laplace transform 
of ip holds true for t — > +0: 

e-*>(A)rfA~^p i r(l-a J )^- 1 + (logt)^g j7 i-i— 1^+ ^ r^+O^" 1 ) 

■°° < j - 1 OsCfc<M-l 

(2.4) 

some coefficients {r&}. 

The proof can be performed, for example, by explicit computation, checking that each 
term in the asymptotics ()2.3|) gives the desired contribution to (|2.4|) . 

Proof of Lemma \2.1\ 1 . Let 

N{\) = ${n | A n < A}, iVo(A) = t{n | A° < A}. 

The main idea of the proof is to approximate N(X) by N (ijj(\)), where ip is a function with 
the asymptotic expansion (J2.1|) . We construct t/> in terms of its inverse as follows. 
The formal inversion of the expansion ()2.1|) has the form 

oo 

A n ~A°+X;^(A°r w , n^oo, (2.5) 
i=i 

where ^ < 772 < . . . and {sj} C R. Fix some sufficiently large M G (0, 00) \ N; we 
have 

An = A° + J2 + 0((X° n )- M ), n^oc. (2.6) 

r/j<M 

Let G C°°(R) be such that 

(i) 0(A) ^ for all A G R and 0(A) = for all A < 1; 

(ii) 0(A) is strictly increasing for A > 1; 

(iii) 0(A) = A + Yl s jA~ % f° r a U sufficiently large A > 0. 

T)j<M 

Let t/) G C°°(0, 00) be such that 0(^(A)) = A VA > 0. Finally, for A > let us write 
iVo(A) = |A + u;(A), where ct»(A) is a 2-periodic function. 
With this notation we have: 



x POO POO 

Y^e~ tx - = e- tx dN(\)=t e- tx N(\)d\ 
1 J —00 J —00 



n=l 



tl e- tx (N(\) - N (i)(\))) d\ + \t I e~ tx i}(\)d\ + t I e~ tx u{i;{X))d\ 



2 

00 * 



F 1 (t) + F 2 {t) + F 3 (t). 



6 



Below we consider separately the integrals F\(t), F 2 (t) and F 3 (t). 

2. Consider F 2 (t). By the construction of ip, we have the asymptotics 

■0(A) = A + Pj X ~° 3 + E q i X ~ j + °( A_M ) A 00 

a.j<M j<M 

with the same exponents and coefficients as in (|2.1|) . By Lemma \2. 21 we obtain: 

m) ~ Yt+ E |r(l-%)^ + i(logt) E *7^y^+ E ^+°(^), * - +0. 

aj<M j<M V >' 0<k<M 

(2.7) 

3. Consider Fi(t). By the construction of 0, we have A n = 0(A°) + 0{n AI ), n — > 00, 
and so 



• ')■ 

n=l n=l 

It follows that F\(t) has at least [M] — 1 continuous derivatives in t on [0, 00) and 
therefore, by the Taylor formula, 

F 1 (t)= E Fi k) (0)t k + o(t [M] ~ 1 ), (2.8) 

0<fc<[Af]-l 

4. Let us prove that F3 has continuous derivatives in t £ [0, 00) of any order, and so 

00 

F 3 (t)~E F 3 (fc) (0)* fc , «-+0. (2.9) 



fc=0 

Fix N £ N. Integrating by parts N times, we obtain 



00 POO 



F ?J {t)=tj e- tx u{<4){\))d\ = t J e-^V(//)u;(//)d// 

./o 

where cu^v is a periodic function. Using the property (iii) of 0, we obtain 

N 

( e -t*(M))(*+i) = e -^){(-t) iV+1 (0 / (/i)) JV+1 +E* Z °( AZ ~ iV ~ 2 ~ 7 ' 1 )}- 



z=o 



It follows that _F 3 has at least AT — 1 continuous derivatives on [0, 00). As N £ N can be 
taken arbitrary large, this proves the statement. 

5. Combining (j22j) - (JUj), and using the fact that M can be taken arbitrary large, we 
get the desired statement. ■ 
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3 Proof of the asymptotic expansion ( 11.31 ) 

1. Asymptotic expansion (jl.lj) . First let us prove that for any bounded from below 
function v G C°°(R), the asymptotic expansion (jl.lj) . (jl.2J) holds true locally uniformly in 
x G R. The expansion (jl.lj) as such is well known, but all treatments of this expansion in 
the literature that we are aware of, assume boundedness of v, whereas here we have to deal 
with potentials of the type x 2 + v (x). Below is a simple argument which shows that the 
boundedness from above condition can be lifted. Let us fix any R > and prove that (jl.lj) 
holds true uniformly in x G [-R, R]. 

Let v G Co°(R) be such that v(x) = v(x) for all \x\ ^ 4R, and let h = —£s + v. The 
heat kernel expansion for C£°-potentials is certainly well known (see e.g. [Zj and references 
to earlier work therein), and so we have 



e- th (x,x) ~ 

V4nt . „ 



^ oo 

= ^2t j a j [v(x)}, 



uniformly in x G [-R, R] ■ Thus, it suffices to prove the estimate 

sup\e- t ~ h ( y x,x)-e- th (x,x)\=0(e- c/t ), t -> +0, c> 0. (3.1) 

\x\^R 

Let xr De the characteristic function of (—R,R) in R, and let <p G C 3O (R) be such that 
4>{x) = 1 for \x\ ^ 2R and <p(x) = for \x\ ^ 3i?. Denoting D = we obtain 



Jo 

ds XRe~ {t - s)h <f>"e- sh XR + f ds Xtfr^fy 'De^XR 

Jo 

[ t ds X Re- {t - s)h re- s ' h XR+ [ dsxRe^-^'De^xn 
Jo Jo 

- [ ds [ d Sl XRe~ {t ~ s)h (j)'De 
Jo Jo 



(^) D2 ve- sJl x R . 



From here, using the explicit formula for the heat kernel e and the well known estimate 
|e- 4/t (x,y)|^-i=exp(-^-tinft;), x,y eR, (3.2) 

we obtain (j3.1j) . 

2. Asymptotic expansion (jl.3j) . Now we are ready to prove the asymptotic expansion 
(jl.3j) . Let R > be sufficiently large so that supp q C (— -R, -R). Let be the characteristic 
function of (-2.R, 2R) and let X2R = 1 — X2R- By the previous step of the proof, it suffices 
to prove that 

Tr( X 2 R (e- tH -e- tHo )x 2R )=0(e-^), t - +0, c> 0. (3.3) 



S 



By (13. 2|) . we obtain 

\\X2 R e~ tH XR\\s 2 = 0(e- c/t ), t -> +0, c> 0; 
||Xi?e-^X2i?lk = 0(e- c/t ), t - +0, c> 0. 
From these estimates and the formula 

X2 R (e- tH - e- tHo )x2R = - f X2 R e'^ s)H XRqXRe- sHo X2 R ds 

Jo 

we get the required result (|3.3j) . 

4 Proof of Theorem 11.31 

We follow the arguments of [3]. First let us assume that A n 7^ for all n. Fix k G N and 
consider formula (Jl.lUj) . The second term in the r.h.s. is meromorphic in C with possible 
poles at s — 1 — |, j = 0, 1, 2, . . . , 2k + 2. The first term in the r.h.s. of (jl.K)j) admits 
analytic continuation into the half-plane Res > —k — ~. As A; can be taken arbitrary large, 
it follows that Z admits a meromorphic continuation into the whole complex plane, all 
poles of Z are simple and located at the points s = 1 — j, j = 0,1,2, ... . 
Next, from the formula 

j /"7+ioo 



we get 



j / 1 7+200 

e" tA " = — / (t\ n )- s T(s)ds, 7 > 0, (t\ n ) E R \ {0}, 
J2^ tXn - ^ = V~ / ~ Ms))t- s T(s)ds, t>0, 7 > 0. 

n=l J-y-ioo 



By a standard argument involving shifting the contour of integration to the left, the last 
formula yields the following asymptotic expansion as t — > +0: 



CO 



5>~ tA " - e~ tX °) ~ E Res ((Z( S )-Z ( S ))rT( S )), t - +0. (4.1) 

First note that Z(s) — Z (s) does not have poles at any of the points s — 0, —1, —2, 
Indeed, if Z(s) — Z (s) did have a pole at s — —n say, then (Z(s) — Z (s))T(s) would have 
a double pole there and then the expansion (J4.1j) would involve a term Ct n logt. But by 
(|1.3jh no logarithmic terms actually occur in the asymptotic expansion. 

Next, by (|1.3jl . there are no integer powers of £ in the asymptotic expansion, which by 
the same argument leads to the conclusion that Z(—k) — Zo(—k) = for all k — 0, 1, 2, ... . 

Finally, consider the case when one of the eigenvalues of H vanishes: A m = 0. Then the 
preceding arguments should be repeated for the sequence {A n }, n G N \ {m}. This leads 
to the same set of results, apart from the formula Z(0) = 0; this should be replaced by 
Z(0) = -1. ■ 
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5 Proof of Theorem ll.ll fi) 

Let us define two solutions = i>±(x, A) of the equation —ip" + x 2 ip = \ip by 

iP° + (x,\) = U(-~,xV2), V-(x,A) = U(-~,-xV2), 

where U is the parabolic cylinder function (see jU §19.3]). For any x G R, the solutions 
il>±(x, A) are entire functions of A. For any A G C, the solutions ?p±(x, A) have the asymp- 
totics 

^° {x, A) = ^°(-x, A) = (xV2)^-^ 2 e- x2/2 (l + o(l)), a; -> +00, 
and the Wronskian u>o(A) = , = (i } -)' x i } + ~ ^-(i^+Yx * s gi ven by 

^) = f ^ = -|r(l±A)cos(f). (5.1) 

At the eigenvalues A° = 2n — 1, the Wronskian wq(X) vanishes and we have 

1>\{x, A°) = (-l)" + V-(^, A°) = 2-^ 2 e- x2 / 2 H n ^(x), (5.2) 

where if n is the n'th Hermite polynomial. 

Next, let ip± = ip±(x,X) be the solutions of the equation —ip" + (x 2 + q(x))ip = A^, 
normalised by 

if)+(x, A) = ^ (x, A), x>supsuppg, 
i/j-(x, A) = ij)2.(x, A), x<infsuppg. 

The eigenvalues A n coincide with the zeros of the Wronskian io(A) = In 
Section |U] we prove the following Lemma, which describes the asymptotics of w(X) as 
Re A — > +00. Let Q be the half-strip 

tt = {A G C I ReA ^ 0, |ImA| 1}, 

for A G f2 let us denote by y/X the principal branch of the square root, so that Re \f\ ^ 0. 

Lemma 5.1. The Wronskian w(X) is analytic in A G Q. The following asymptotic expan- 
sion holds true: 

ro(A) . _| r(¥) ( cos(f ) £ + M f ) g ^) , (5.3) 

as I A| — > 00 , A e O. i^ere Qj, Pj G C are some coefficients, Qo = 1, Pq = 0. 

Given Lemma 15.11 we can prove Theorem II. If i) as follows. Fix any sufficiently small 
e > 0, denote B n>£ = {z \ \z — A°| ^ e}, and let r n)£ be the contour dB ny£ oriented anti- 
clockwise. By Rouche's Theorem combined with a simple continuity argument, we obtain 
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that A n G B n £ for all sufficiently large n. Next, the zeros of w in the half-strip Q coincide 
with the zeros of 

y/nw(X) 



w(X) 



2r(i±A)" 



It follows that for all sufficiently large n we have 

K = ±[ X^dA. (5.4) 
2m J Tn e w{X) 

By analyticity of w, the asymptotic expansion (|5.Hj) can be differentiated. Thus, we obtain 
the following asymptotic expansion for A G r n e , n — > oo: 

= --tan T + — ^(A) + ^^ 1 (A)tan T + ^(tan T r(v / A)- m+ V(A), (5.5) 

where the functions g m (X) are analytic in A G VL and have the asymptotic expansions 

oo 

9mW~'52c mk (V\)- k , |A| -> oo, A efi. (5.6) 

fc=0 

Substituting the expansions ()5.5|) and (|5.6|) into ()5.4|) and computing the integrals of the 
type f T (tan 2 y-)" 1 A J <iA, we arrive at the expansion (jl.4j) . ■ 

6 Proof of Lemma 15.11 

Let x > sup supp q; then 

w(X) = (Vi_(x, A))^° (x, A) - ^_(x, A)«(x, A));. (6.1) 

We will use this formula and construct ip- in a standard way as a solution to the integral 
equation 

ij.(x,X) = ^°_(x,X)+ [ G x (x,y)q(y)^(y,X)dy, (6.2) 



where the integral kernel G\(x,y) is given by 

G x (x, y) = ^r«(x, X)^°_(y, A) - ^ _{x, A)<(y, A)). (6.3) 

100(A) 

The kernel G\(x,y) is an entire function of A due to the analyticity of ip±(x,\) and the 
relation (I5~2l . 

Let i? > be sufficiently large so that suppg C (—R,R). Denote A = [— 2R, 2R]; let 
L x : C(A) — ► C(A) be the Volterra type integral operator from (|6.2j) . 

/x 
G\(x,y)q(y)f(y)dy. 
2R 
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Then the solution of the integral equation (J6.2j) can be written as 

oo 

n=0 

and so for the Wronskian (|6.1|) we have the series representation 

oo 



IV [ 

n=0 



Lemma 6.1. For any n G N and any x G (R,2R), the Wronskian W(L^ip^_,ij}^.)(x) is 
analytic in A G Q and one has the estimate 

\W(Lir-,r + )(x)\^^^\TC-t)\, C(A) = 0(|A|- 1 / 2 ), |A|^oo,AGfi. (6.4) 
The asymptotic expansion 



oo q{u) oo p(n) 



W(Ll^° + )(x) ~ r(*±*) ( cos(f ) £ + sin(f ) g ) , (6.5) 



with some coefficients Q ™ , Pj holds true as |A| — > oo ; A G O. 
Clearly, Lemma 15.11 follows from Lemma 16.11 

Proof of Lemma \6. 1\ 1. It is convenient to introduce two linear combinations e+ and 
e_ of of the solutions ip±'- 

/=o(l-A)/4 

A) = 7" + (e-^y^ + (x, A) + e^'Wi A)) , 

e_(x, A) = e+(— x,X). The solutions e±(x, A) are analytic in A G f2 (with removeable 
singularities at A° — see (|5.2Jl ). These solutions are chosen so that they satisfy the following 
asymptotic expansions: 

e ± (x,A)~e ±i ^'^l + ^-^^J, A^oo, A G Q, (6.6) 

(oo J^-(x)\ 
±iV\x + Y,jj=^J, A^oo, AGfi, (6.7) 

where Rj, Rj are polynomials in x. The expansion (|6.6|) follows directly from the formulae 
19.9.4, 19.9.5, 19.4.2 of pQ, and (J6.7|) is obtained by application of the recurrence formulas 
P §19.6]. 
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2. Let us first prove the bound (|6.4j) . We have 

W(LW°_^l)(x) = (L^°_(x,X))' x ^ + (x,X) -LW _(x,\)(i> + (x,\))' x ; (6.8) 

let us obtain appropriate bounds for each term in the r.h.s. of (16. 8|) . Expressing in 
terms of e±, 

4>l{x, A) = ^- 2 ^ 1 )/ 4 r(l±A) (^V, A) - e-^ +1 )/ 4 e ± (x, A)) , (6.9) 

and using ()6.6|) . ()6.7|) . we obtain 

||^(-,A)|| C (A) = 0(|2 A / 4 r(i±A)|), |A|-oo, Aefi, (6.10) 
ll(^(-,A));|| C ; (A) = 0(|A 1 / 2 2^r(i±A)|), |A|-oo, AgO. (6.11) 

Next, expressing the kernel G\(x,y) in terms of e±, 

1 r(^) 

G\(x, y) = — ^^ x ^ (e + (x, A)e_(y, A) - e_(x, A)e+(j/, A)), (6.12) 
and using the asymptotics ()6.6|) . we obtain 



sup sup 

|a|^R|y|<fl 



G? A (z,y)| = 0(|A|- 1 / 2 ), |A|-oo, Aefi. 



Using this estimate and the fact that L\ is a Volterra type operator, we obtain 

MkAW7(A) ^ C(A) = 0(|A|- 1 / 2 ), |A|-oo, AGO. (6.13) 



n! 



Finally, in order to estimate the term (L^i/;^_ )' x7 let us introduce the operator L' A : C(A) — > 

L' x : /(*) ~ jT dGx £ y) q(y)f(y)dy. 

Then (Z#0°(ar, = £'a l a~V-(^ A). Using the asymptotics IfKlfl). (l6~7j) . we obtain 

\\L' x \\ C (A)-+c(A) = 0(l), |A|->oo, AGO. (6.14) 

Combining (jFTSjl . (l6~TUj) - (|6~T4T) . we obtain (IQll . 

3. Let us prove the asymptotic expansion ()6.5|) . Using ()6.9|) . we obtain 

, O = i- 2 ( A - 1 )/ 2 r(i± A ) 2 (-lU(^e + , e + ) - lU(L^e_, e_) 

+ e -^ (1+A)/2 lU(L^e_, e+) + e in{1+X)/2 W(L n x e + , e_)) 

Denote 

^(x,A) = L "^ ) ; (6.15) 
e±(x, A) 
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by 1)6. 6 Jl . the denominator does not vanish for all sufficiently large A. Using this notation, 
we obtain 



+ (^ + (^A));o(|r(i^)|) + (^(x ) A)y x .o(|r(i±A)|). 



It suffices to show that have the asymptotic expansions 

Sf(x 

(Va> 



00 S ± (x) 

^(x,A)~£-^4, |A|-oo, XeQ (6.16) 



j=n 



for some coefficients £?■ G C(R), and that for any x G (R, 2R), 

(g±(x,\))' x = 0(\\\-°°), |A|-oo, Aefi. (6.17) 
By the definition of g„, we have 

9n+ii x A) = [ -^r 1 tt gn(y, A)e ± (y, X)q(y)dy. (6.18) 
i-i?e±(x,A) 

Using this formula, the expression (J6.12)) for G\(x, y) and the asymptotics ()6.6j) . the expan- 
sion ()6.16j) can be easily proven by induction. The asymptotics (16.17)1 follows by differen- 
tiation of (I6~T51) . 
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